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^ ' Abstract 

^ . We consider Schur function expansion for the partition function of the model of normal 

matrices. We show that this expansion coincides with Takasaki expansion [Sj for tau 
' functions of Toda lattice hierarchy. We show that the partition function of the model of 

normal matrices is, at the same time, a partition function of certain discrete models, which 
can be solved by the method of orthogonal polynomials. We obtain discrete versions of 
— ' various known matrix models: models of non-negative matrices, unitary matrices, normal 

matrices. 



> 



o 

)Q 1 Introduction 

o . 

' The model of normal matrices (MNM) was introduced and studied in pp, j2] and [3]. A 

I • matrix M is called normal, if it commutes with hermitian conjugated matrix: 
^ ; 

=3 ■ [M,Mt]=0 (1.1) 

>- [ One can bring the matrix M to its diagonal form via the unitary matrix U: M = UZW , where 



X 



Z = diag(zi, . . . , Zn), with the eigenvalues of M. It is clear that ikfl" = UZU\ where the 



i-i diagonal matrix Z is the complex conjugate of Z. 
— . The model of the normal matrices is defined by its partition function as follows: 

Zn = C J cin(M)eT'-^^(*^)+T^^^(*^')+T'-^(^'*^') , (1.2) 

where C is a normalization constant. The integration measure is defined as 

n 



J2 



i=l 



where d^:U is the Haar measure on the unitary group U{n), and A{z) is the Vandermonde 
determinant: 

n 

A{z) = det{zr%,,=,_^ = - z,) 

i<k 

For the case n = 0, one sets A(z) = 1. 
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The potentials Vi{M), V2(M"I') are defined by their Taylor series as follows: 



m=l m=l 

where t = {ti,t2, ■ ■ ■), t' = (t[,t'2, . . .); tm,t'^ are called coupling constants. The interaction term 
V, a function of two variables, is not specified here. Integral ()1.2|) is supposed to be convergent; 
in the present paper, we consider it as its perturbation series in the coupling constants. 
After the change of variables M (f/, Z) and the integration over U{n), one obtains 

Z^ = Z^{V;t,t')= [ ■■■ [ |A(z)|2TTe^i(^^)+^^(^">)e^(^-^"»)rf2^, , (1.3) 
Jc Jc -^^ 

where the integration is over the complex planes of eigenvalues Zi. Here we choose the normal- 
ization constant C = C{n) in such a way that C ■ Vol(f/(?2)) = 1. 

The model of normal matrices has applications in a description of quantum Hall droplets. 
As it was shown in the particular case of the model, namely, the case when V{z, z) = —\z\^ 
in its n — s> oo limit, is related to the interface dynamics of water spot with a constant source 
surrounded by oil; in addition, it is related to certain old problems of complex analysis. The 
relation of these problems to Toda lattice equations [T3] and its dispersionless analogue [Sj was 
found and discussed in jlj. 

The integrable structure of the model of normal matrices was found and discussed in jH]. 
The coupling constants t, t' and the size of matrices n were identified with the so-called Toda 
lattice higher times. 

Remark 1 One can absorb both Vi and V2 into the term V . We prefer not to do so, since 
we are interested in the integrable structure of the model; we therefore keep the parameters 
t = (ti,t2, ■ ■ ■), t' = (t[,t2, . . .) as the Toda lattice higher times, see below. 

Remark 2 // the term V enjoys axial symmetry 

V{z,z) = V^,i\z\') , (1.4) 

then, for a certain choice of coupling constants, there exists a different integrable structure: 
partition function / li.jj)) is an infinite- soliton tau function of a "dual" TL hierarchy, higher 
times of the dual hierarchy being related to the moments of e^™ . See fT^ about this and related 
topics. 

Schur functions and partitions Polynomial functions in many variables, like the Schur 
functions, are parameterized by partitions. A partition is a sequence of non-negative integers 
in the non-increasing order: 

A = (Ai,A2,...,A„...) , Ai > A2 > ... > A, > (1.5) 

where we identify (Ai, . . . , A„) with (Ai, . . . , A„, 0), and in the case of infinite sequence, assume 
Ar = as r ^ 0. The Aj in ()1.5|) are called the parts of the A. The number of nonzero parts 
is the length of A, denoted by i{X). The sum of the parts is the weight of A, denoted by |A|. If 
n = |A|, we say that the A is a partition of n. The partition of zero (where Ai = 0) is denoted 
by boldface 0, to distinguish it from number 0. The set of all partitions, including 0, is denoted 
by P. 



The (Young) diagram of a partition A is defined as the set of points (or nodes) G Z^, 
such that 1 < J < Aj. Thus, it is a subset of a rectangular array with i{X) rows and Ai columns. 
We denote the diagram of A by the same symbol A. For example, 



;i.6) 



is the diagram of (3, 3, 1). The weight of this partition is 7, the length is equal to 3. 

The partition whose diagram is obtained by the transposition of the diagram A with respect 
to the main diagonal is called the conjugated partition and denoted by A*. 

The product of hook lengths H\ is defined as 

= n ^i^' = Ai - i + A* - J + 1 , (1.7) 

where the product ranges over all nodes of the diagram of the partition A. 
Given number g, the so-called hook polynomial Hx{q) is defined as: 

Hx{q) = n (1 - ^'")' = A, - z + A* - J + 1 (1.8) 

i,jex 

In what follows, we also need notations: 

k 

n(A) := 5^(2 - l)Ai , (1.9) 

i=l 

P _|_ 777,^ 

(a)A := (a)Ai (a - 1)a2 ■■■(«- ^ + l)Afc , {a)m ■= — — , (1-10) 

(g^ g)A := {q^; q)xAq^-'; q)x, " " " {q"''^': ?)a. , (g'^; q)m := (1 - g'^) ■ ■ • (1 - g'^+"-^) , (1.11) 

where k = i{X). We set (a)o = 1 and (g"; g)o = 1. 

We now consider a semi-infinite set of variables t = (ti,t2,^3, • • •)• Given partition A, the 
Schur function sx(t) is defined by 

oo oo 

SA(t) = det(/;,A,-i+i(t))i<j ,,<^(A) , where ^ z*^/ifc(t) = exp ^ z^t^ , (1.12) 

fe=0 m=l 

and, for < 0, we put = . The hk{t) is called the elementary Schur function. 

There is another definition of the Schur function; it is the following symmetric function in 
the different variables ), where n > i{X): 

= ''^fZ?''"'" (1-13) 



for the zero partition one puts Sq{x) = 1. If 

n 

t = t(x(")) = (ti(x(")),t2(x(")), . . .), t™(x(")) = - 

m < ^ 



m 

1=1 



then definitions ()1.12|) and ()1.13p are equivalent jH]: 

SA(t(x(")))=.A(x(")). (1.14) 



Remark 3 From definition M.l'^) it follows that SA(t(x("))) = if > n. 



The Schur functions Sx{xi, . . . where £(A) < n, form a basis in the space of symmetric 
functions in n variables. We use the underhne in only to distinguish the two definitions. If 
an n X n matrix X has eigenvalues xi, . . . , we may denote s^^{xi, . . . , Xn) by sx{X), without 
underline, since in this paper the Schur function with uppercase argument is used only in this 
sense. 

Remark 4 If any of Xi coincides with any ofyi, say, x^ = yj, then we assume that the function 
qV{z,z) yQfiigJiQg at z = x~j^^ so that the integral M.'^) is convergent. 

Soliton theory While the KP hierarchy of equations, jJUlj is the most popular 

example in the soliton theory, the two-dimensional Toda lattice (TL) is another important 
equation, first integrated in P^, and carefully studied in in the framework of [T^ : 

dtA'^i^n = e'^"+i-'^" - e'^"~'^"-i (1.15) 

This equation gives rise to the TL hierarchy which contains derivatives with respect to the 
higher times ti, t2, ■ ■ ■ and t[, t'2, . . . . This equation was applied to the study of models of 
random Hermitian matrices (the so-called one- and two-matrix models) in JH]; 

The key point of the soliton theory is the notion of tau function, introduced by Sato (for 
KP tau function; see [T2]). The tau function, a sort of a potential which gives rise to the TL 
hierarchy, is a function of the two sequences of higher times, ti,t2, • • • and ^'1,^2, . . . , and the 
discrete variable n: r = r(n, t,t'). More explicitly, we have ^21, [T^ : 

u = 2dl In r(n, t, t') , ip^{t, t') = In ""^"^ + 



r(n,t,tO 

Power series solutions The formal power series solutions of the TL hierarchy were described 
by Takasaki in the form of double series in the Schur functions over partitions [Hj (see also P , 

r(n,t,t')= ^ (7A,A'(^)sA(t)sv(t') , (1.16) 
A,A'eP 

where the coefficients ^a,a' can be presented as certain determinants or, alternatively, just solve 
special bilinear equations ^2]- We are interested in the semi-infinite Toda lattice, i.e., when 
the index n runs over the set of nonnegative integers Z>o. In the semi- infinite case, we have 
the restriction 

£(A),£(A') < n 

of the summation range on the right hand side of p.lfij) . By Theorem 1 of in this case 

g x,x' (n) = det [Gh^h') 

\ ■' / i,j=l,...,n 

hi = Xi-i + n>0 , h'i = X'i-i + n>0 , 

and semi-infinite matrix G is related to the factorization problem \E'^(t, t') = \E'~(t, t')^, where 
the semi-infinite matrices \E'^ are Baker- Akhiezer functions [15j. 

If all but finitely many terms in (jl.l6j) vanish, it gives a rational solution, in which case the 
formula may be best known. We are mainly interested in the case where fjl.l6p gives an infinite 
sum. 



Tau functions of hypergeometric type In the sequel, we shall also consider the diagonal 
case: g\\'{n) = if A 7^ A'. In this case, generically (i.e., if go^o{n) 7^ (Vn); this restriction is 
not essential) we have g\\i{n) = 6x,x'rx{n)go^o{n), where the coefficients rx{n) are given in the 
following way: 

Consider a function r = r{n) of a single variable n E 1^. Given partition A, we define 

rx{x) = J]^ r(x + J - i) , ro = 1 . 

Namely, we place r{x+j — i) at each node of diagram A, and multiply them over all nodes. 
The "content" j — i oi node («, j) vanishes on the main diagonal. For instance, if A = (3, 3, 1), 
then (see the diagram in p.6|) ) rx{x) = r{x + 2){r{x + l))^(r(x))^r(a; — l)r{x — 2). 
We define the tau function of hypergeometric type as 

r,(n, t,t') = J2^x{n)sx{t)sx{t') , (1.17) 

A 

which solves the TL equation ()1.15|) . rewritten in the form 

dtA[^n = r(n)e'^"-^-^" - r(n + l)e'^"-'^"+i , = -0, + e„ , (1.18) 
where we set r{k) = e^'=~^'=-i and 

= r.(n,t,tO 

When r(0) = 0, we obtain semi-infinite Toda lattice ()1.18p . In what follows, when we consider 
semi-infinite Toda tau function of the form t^, we always assume r(0) = 0. 

A simplest example yields a useful and simple relation (bosonic (left hand side) and fermionic 
(right hand side) representations of the vacuum TL tau: 

00 

exp ^ m-im-i'm = ^a(7)sa(7 ) , (1-19) 

m=l AeP 

which is a generalized version of the Cauchy-Littlewood identity jO]. Here we use notations 
7 = (7i)72, • • •) and 7' = (7^,72, ■ ■ ■) in place of t and t', for future purposes. (Since we did 
not find this identity in literature, we presented a proof in (23j). 



2 Discrete matrix models 

By discrete matrix model we mean models of random matrices, with eigenvalues lying on a 
lattice, the type of the lattice is different for different models and will be specified later in 
the section "Associated discrete models and discrete one-matrix models". We shall consider 
discrete counterparts of the following models, where integrals over eigenvalues are replaced by 
sums. Namely, we shall consider the following replacements. 
Model of normal matrices: 

n n 

\A{z)\^ n e''^(^')+''^(^"^)e^(^-^">)d2z, ^ ■ ■ ■ \A{z)\' J] ^v,(zO+v,(-z.)^v(z.,2.) = z, (2.1) 

1=1 i=l 



Model of two positive Hermitian matrices: 

„ n n 

/ A{x)A{y) Yl e^^("^)+^^(^')e^("-^')(ixi% ^Y.'"Y1 ^(^)^(?/) H e^^^^'^+^^^^^^e^^"-^') = Z2 

•^^+' i=l i=l 

(2.2) 

Let us notice that the well-known standard two-matrix model is described by specific choice of 
the interaction term: V{x,y) = xy. 

Model of positive Hermitian matrices: 

P n n 

/ A(x)2j]e^(^»)d2;,^ J]--- J]A(a;)2j]e^(^') =^3 (2.3) 

•^^+ i=i i=i 

Model of two unitary matrices 

A(x)A(y) JJe^^^^^^+^^'^^'^e^^^'*)^^ ^ ■ • • A(x)A(y) J] e'^^^^^^+^^^^^^e^^^-^') = Z, 

i=i ^^y^ t=i 

(2.4) 

where, for the standard case, one takes V{x^y) = x ^y ^. 
Model of unitary matrices 

/n n 
|A(^)|2 J]e^^(^»)+^^(^"»)e^(^-^"»)rf2^, ^ ^■•■5^|A(z)|2 J]e^^(^»)+^^(^"»)e^(^-^"») = Z, (2.5) 
i=l 1=1 

Kontsevich-like model: 

/ A{x)]le^^y^e^^''^Uxi^y2---y,^(^)I\^'^'''^^^'''^ = ^<^ (2-^) 

The multiply integrals in the left hand sides are integrals over the eigenvalues of the related 
matrix models. The cases of the generalized interaction term V{x, y) in ()2.2p and ()2.4p should 
be considered separetely in order to be considered as integrals over matices. We do it in section 
"Two- matrix models with generalized interaction term" below. 

What we are going to prove is that by choosing parameters t, t' in a special way we obtain 
ZniV; t, t') = patition function of discrete matrix models , Zi, . . . , Zq (2.7) 

3 Schur function expansion of the partition function 

In [22], for the Gaussian interaction term V{z,z) = — the following perturbation series 
was obtained 

Z„(V;t,t') =n! HASA(t)sA(t') , (3.1) 

AeP 

£(A)<n 

where {n)x is defined in At the present paper we generalize this result. 

Given V{z, z), let us introduce bi-moments gkm, where /c, m = 0, 1, 2, . . . , as 



dkr, 



^fc^mgF(.,f)^2^ , (3.2) 



assuming the integral on the right hand side is well-defined. 



Theorem 1 The partition function of the model of the normal matrices il.!^} has the following 
perturbation series in the coupling constants t and t' : 

gx,\'{n)sx{t)sx'{t') , (3.3) 

A.A'eP 
e{\),e(\')<n 

where 

g \,x' (n) = det (gh,h') , (3.4) 

hi = \i — i + n>0, h[ = \'i — i + n >0 (3.5) 

Proof. The proof repeats one of the proofs of ()3.1|1 . written down in [22j. Taking 7^ = tm, 
7^ = (1/m) TrM™", and applying to the term e"^^^ , and applying it again to the term 

e^^2^ with 7„ = t'^, = (1/m) Tr(Mt)'^, we obtain 

Zn{V -t.t') 

„ n 

= E sxit)syit') / det (..'^0 det (z/^) He^^-^-^^d^ , 

A.A'eP "^^ - ' - 1=1 

^(A),^(A')<n 

where we used definition ()1.13p to get the determinants in the integrands; these determinants 
are parameterized by the sets {hi \ i = 1, . . . , n} and {h'^ \ i = 1, . . . ,n}. The rest of the proof 
is quite similar to the calculation suggested by John Harnad for |22j . 

Given a pair of partitions (or, by ()3.5p . given a pair of sets {hi,}, {h[}), we develop each 
determinant as a sum of n\ monomial terms: the first determinant yields the sum of monomials 

over all permutations a of numbers l,...,n: Ylaes i^T^^^^"^^ YYi=i ^i'''''^ each monomial is 
parameterized by the set of non- negative numbers {hi} and by an element of the permutation 

h' 

group a. Similarly, the second determinant yields '^cri—Y^^^^'^^Yli^i"''^''^ each monomial is 
parameterized by the set {h[} and by a. After integration, each pair of monomials produces 
the product of n integrals of type ()3.2j) . Notice the restriction ^(A), ^(A') < n, which follows 
from Remark 121 Given a pair of sets {hi} and {h'^}, we gather all the terms into determinant 
fl3.4|) and obtain the Theorem. 

Note that the very first term of (|3.3p yields the formula presented in [S] 



Z„{V; 0, 0) = n\go,o = n! det (^j z''-'z''-'e^'^'''^d^z^ 



i,j=l,...,n 

Corollary 1 Using Theorem 1 of we get another proof of the fact that the partition function 
of the model of normal matrices is a tau function of the (semi-infinite) TL hierarchy, 

Inverse moment problem One can ask what kind of potential V{z, z) of MNM p.2|) pro- 
vides a given form of potentials in the discrete matrix models we have obtained. We need to 
present the following comment. 

Remark 5 Since the moments gmn (^re the higher derivatives at the origin of the Fourier 
transform of , and since every formal power series is the (formal) Taylor series at the origin 
of some C°° function, one can, in principle, regard the moments gmn as independent variables. 
Having this fact in mind, in the future we may write each bi-moment gmn of \3. 2) in a form 
analogous to e^'^'^'^^e^ in il.!^} : 

gum = e^^('=)+^^("^)e^- , (3.6) 
and assume certain kind of parameter dependence for Vi{k) and V2{m). 



Note that, since the Taylor series of a function does not determine the function uniquely, 
the moments {gmn}m,n>o do not determine uniquely. Moreover, in popular function spaces 
like or Sobolev, the set of all e^'s whose bi-moments coincide with given {gmn} is not closed 
and may even be dense, so it is hard to give an explicit method or formula to reconstruct one. 

Nevertheless, among those there is at most one which decays exponentially, i.e., 0{e~^^^^) 
as \z\ — > oo, since the Fourier transform of exponentially decreasing function is real analytic, 
and hence is determined uniquely by its Taylor series. 

Parameter dependence of {gkm}, like in ()3.6|) . may often be introduced easily: Suppose, as 
a function of r G M and w = e^^ E {\w\ = 1}, e^^^'^-* = e^*^^"''™ has analytic continuation 
in the ly-variable to C \ {0}, and the resulting function still is rapidly decreasing in r for 
any w G C \ {0}. (This is slightly more general than the axial symmetric case. It includes, 
e.g., polynomials in z and z with coefficients being rapidly decreasing functions in \z\.) Let 
gVfcm ._ j^^k^m^v(z,z)^2^^ Then for any gi, g2 > we have 



hence, introducing two sets of parameters t = (ti, ^2, • • •) and t' = {t[, t^, ■■■)■, and setting 

i,j>0 



where D := —zd/dz — 1 and D := —zd/dz — 1, at least as a formal power series in t and t' we 
have 



A simple diagram summarizes what we get: 



^V{z,z) deformation^ 



bimoments 



bimoments 



(3.8) 



^V^^ deformation^ e?(t,g?)+5(t',gr)e^'=- 



It is more illuminating to put e^^^''^?'''^^^'^' ''^^^ e^^^'^^ in place of e^*-*'*-^'^-* in this diagram, but the 
construction ()3.7p may not work if is in more general form, or if the condition gi, g2 > is 
not satisfied. Nevertheless, we can rely on less explicit argument in Remark El to assume the 
existence of e^'>*' subject to ()3.8p . 



Axial-symmetric interaction term If we take axial-symmetric interaction term ()1.4|) . 
V{z,z) = Vaxi(|-2p), then gkm = Skmgmm is a diagonal matrix, which we parametrize by a 
set of numbers ^i; ^2, • • • as follows: 

r+oo 

gmm = Tr / a;™e^-'(^)da; = e^™-«» , (3.9) 
Jo 

so that 

n 

1=1 

If both A, A' are zero partitions, we have (70,0 = Y[n7=o9mm- Then we obtain 

^n(Kxi; t, t') = go,on\ J2 e^'=^ «"+^-^-«"-^SA(t)sA(t') (3.10) 



Let us notice that a number of matrix integrals (including the generalization of Itzykson- 
Zuber, Gross- Witten and Kontsevich integrals, and various integrals over complex and rectangle 
matrices) may be written in form of series ()3.1U|) . where the variables ^ are specialized via the 
choice of matrix model, see pO] . 

In the following sections, we consider the right hand side of (jH.Hj) and of as models 

Zi, . . . , Ze, where sums over numbers h = {hi, . . . , hn) and h' = {h[, . . . , h'^) play the role of 
sums over eigenvalues. 

4 Specialization of coupling constants. Associated dis- 
crete models and discrete one-matrix models 

Using the Schur function development ()3.3|) . we shall show that there exist associated discrete 
models of random matrices, whose partition functions coincide with the partition function of 
the model of the normal matrices. 

We introduce the following notations: 

too = (1,0,0,0,...) , (4.1) 

/a a a \ , , 

t(a,l)^ {-,-.-...), (4,2) 

t(oo,g) = (ti(oo,g),t2(oo,g), . . .), t^(oo, g) = — — , m = l,2, (4.3) 

m[i q ) 

1 - (g")™ 

t(a,g) = (ti(a,g),t2(a,g),...) , t„,{a,q) = — — , m = l,2, ... (4.4) 

m[l — g™) 

These choice of times was also used in |18j |21j. 

Note that t(a,g) tends to t{oo,q) (resp. t(a, 1)) as a ^ cxd (resp. q 1). As for too, if / 
satisfies /(cti, CH2, cH^, ■ ■ ■) = c'^fiti,t2, ts, . . .) for some G Z, we have ^'^/(t(cx3, q)) /(too) 
as h := In q —>■ 0. 

Lemma 1 For a partition X = (Ai, A2, . . .), let hi := n + Xi — i {I < i < n), where n > i{X)). 
Then 

^A(too) = ^ = , (4.5) 

Hx [li=i hi\ 

gA t a,l = — - = M — — — -T — , 4.6 

Hx Ui=ihilfJi T{a-t + l) 

n n 

A{h) := l[{h, - h,) , A{q') := JJl?"' - , 

i<j i<j 

where for H\, H\{q), n{X), {a)x and {q"';q)x see respectively 7| ), / ti.<^) . / ti..9j) . and 
m.ll]} . Note that those quantities ^.^ - ^.^ are independent of the choice of n > i{X). 



Relations (j4.5p and (j4.6p can be found in [U' (some notations are different) for n = i{X), 
and we only need to add that if we have ()4.5|) and other formulas in the Lemma for the value 
n = i{X), then we have it for all n > i{X), as seen by a straightforward calculation. First 
equalities in the relations ()4.7|1 and ()4.8|1 can be easily obtained from well-known relations of 
[Oj. The second equality in ()4.8j) follows from known formulae of [Oj. 

Using this Lemma, we consider the Schur function development of the partition function of 
the model of the normal matrices ()3.3|) . Choosing the higher times according to ()4.1|) . to ()4.2|) . 
to ()4.3|) . and to ()4.4|) . we obtain the following four models (A), (B), (C) and (D), respectively: 



90,0 

n 



(A) Hermitian discrete matrix models Z2, with positive eigenvalues Consider the 
case 

t = (ti,o,o,...) = titoo, t' = (t;,o,o,...) = t;too (4.9) 

for model ()L2|1 . and use sx{ttoo) = t^^^sx{too)- By Theorem 1 we have 

Z^{V;t,t') = C I rff](M)e^^^(^^'^^^)+*^^^'^+*i^^^^ = n\ ^ t\^\t',r\gxAn)sx{t^)s,it^) 

\,x' 

e{\)/{X')<n 

Note that this integral is a kind of Fourier transform. 

For a pair of partitions A, A', take any n G Z such that n > i{X), i{X'), and let hi := n+Xi—i, 
h[ := n + A'j — 2, for 1 < i < n. Then by using ()4.5j) . the perturbation series for the partition 
function of the model of the normal matrices becomes 

■^ra(^) too, too) 

n 

= ^'-90,0 V V A(/.)A(/^One"'"'^''^'^"''^'^'^^'^det(^?,.,,) 

hi>...>h„>Oh[>...>h'„>0 i=l 

n 

f J2 A(/^)A(/.')ne"'°'^'^^'^'"^'^'^^'^det(^,,,,,)^ _^ ^ , (4.10) 

h,h'>0 i=l 

where we got the factor 1/n! instead of n\, as we switched from the summation over the cones 
hi >■■■> hn > 0, h'l >■■■> h'^ > to the summations over all non-negative integers 
hi, . . . ,hn and h[, . . . , h'^, denoted by X]/i>o Sh'>o throughout the text. According to 

f)3.2j) . the factor det igh.h'. ) is defined by the interaction term V in the model of the 

V ^/ i,j=l,...,n 

normal matrices and describes the interaction between the random matrices of the discrete 
model. 

Using the skew-symmetry of A{h)A{h'), one can replace the determinant on the right hand 
side of ()4.im) by the product of n terms ghih[, i = 1, . . . , getting n\ as a multiplier. 

Recall that |A| = "^^hi — (n^ — n)/2. Finely, we find that the perturbation series for 
the partition function of the model of the normal matrices is equal to the following partition 
function of a discrete model: 

Z„(\/;t,t') = (M;)^^o,o J2 A(/.)A(/^')ne~'°'^'^'''^"'''^'^"''^'''''''*'^'''''*'^^''''''''^ ' 

h,h'>0 i=l 

(4.11) 

where each sum ranges over the n-tuples of all non- negative integers, hi, . . . ,hn and h[, . . . , h'^, 
respectively; the times t,t' in the left hand side are chosen according to ()4.9j) . and we remind 
that 



n It It 

j zH'^e''^'~'U\ , A{h) = l[{hi - hj) , A{h') = lliK - h'j 



We notice that the first non-trivial term in sum (j4.1ip . corresponding to hi = h'- = n — i, is 
equal to one. 

In case gkm = f{km)e'^^^''^^^'^^"^^ (see Remark Ej) the sum ()4.11|) is a discrete analogue of 
the model of two Hermitian random matrices with non-negative eigenvalues with generalized 
interaction term Z2, see ()7.13|) below. In case f{km) = e~^"^ we obtain a discrete analogue of 
standard model of two Hermitian random matrices. 

If we take axial-symmetric interaction term ()3.10p . we obtain the partition function of the 
discrete one-matrix model of non-negative matrices Z^: 

^n(Kxi;too,too) = (M;)"^^70,0$^A(/.)2n^'''°''^'''°^*''^^^^'""^""' (4.12) 

h>0 i=l 

(B) We obtain different versions of models Z2, Z^ if we choose 

t = (xa, x^a, x^a, . . .) , t' = (ya, y'^a , y^a , . . .) (4-13) 
Then, (as done in [20]) one can see that p.2|) takes the form 

Zr,{V;t,t') =C j c/^](M)eT■•^(^^'*^')det(J„-xM)-"det(/„-^/Mt)-"' 

Using results of jH] and the definition of hi, we find that, for t as in ()4.13j) . 

x\M ^Y{a-n + hi + l) 

sa t = — M — — -— — 

H\ fji T[a-i + 1) 

Now, for t,t' of form ()4.13|) . we get a partition function with the potential terms slightly 
different of the similar terms in (14.1111: 



r(a-n+fe^ + l) , , r(a'-n + fe^ + l) 

h,h'>0 i=l 



Z^{y-Xt') = {xyy-^go,o A(/.)A(/^0 H 



For axial-symmetric case ()1.4j) . we obtain the following partition function of the discrete 
model of random non-negative matrices: 

^n(Kxi;t,t') = (xy)"^^o.O^A(/^)^n^ r(.--..)r(.^W -^(^-0 +>^.H^y)Hu^-,.-. 

h>0 i=l 

(4.14) 

where variables are related to the axial-symmetric V via (j3.9p . 

Another example: Take the axial symmetric Gauss interaction term Vaxi = — Tr(MM'''). 
Then 



Z„(Kxi;t(a,l),U = C fdQ{M)e-^^'''''^+'^^''det{Ir.-yM^r^ = go,oJ2(yt 



,|A| 



HA(ci) 



A 



HI 



h>0 k=l 



(C) Discrete models related to unitary, ^4,^5, and to normal matrices Zi Choosing 
the times as in ()4.3|) . we obtain that the partition function of the MNM ()1.2|) takes a form 

<m{M) 



J ' 



det(a;M; gi)oo det(?/Mt; 52)00 
where we use the notation for m = 00, and where 

1 x"^ 1 1 , 1 y"^ 1 1 

t ~ {-^^ ) TT 1 2 ' ~o~ 1 ■ ■ ■) 5 t = [y- 



Using the arguments similar to the previous cases, we obtain the following discrete model 

n 

h,h'>0 i=l 

Let us take qi = q2 = Q and |g| = 1 (then we set q = e^^^). 

In case gkm = f{q''q"^)e^^^'^ )+^2(ij (^ggg Remark E)) the sum ()4.1H) is a discrete analogue 
of the model of two unitary random matrices with generalized interaction term Z4, see fj7.13p 
below. In case obtain a discrete analogue of standard model of two unitary 

random matrices. 

For axial-symmetric interaction term V ()1.4|1 in case qi = q2 = q and using parametrization 
()3.9|) . we obtain the discrete model one- matrix model: 

n 

^n(Kxi;t,t') = (xy)^<7o,o$^ |A(g'^)pn^"''°^'^'^''^'"''^'^'^"'^''"^""' ' (4.15) 

h>0 i=l 

which, in the case where |g| = 1, can be interpreted as model of unitary matrices Z^, and in 
the case where |g| < 1, gives a discrete version of the model of the normal random matrices Zi, 
where instead of the integration over the complex plane, we have summation over the spiral 
{g'^|/i = 0,l,...}. 

(D) This case is rather close to the previous one. The choice of times as 

/ 1 — qi x'^ 1 — ql"- x^ 1 — qf"- \ , / 1 — gg' 1 — ^ ~ ^l"' 

t ~ ( -^"j 5 Tr~i 9"' ~o~~i • • • ) 5 ^ ~ [ y' 



. 2 l-g?' 3 1 - g?' ' ' ' V" 1 - ^2 ' 2 1 - gi ' 3 1 - gi ' ' ' ' 

in the MNM model (jl.2p yields the integral 

{xqtM; gi)ooWMt;g2)oo 



Z4V;t,t')=C J ^^^](M)e^^^(^■^^'Met 



(xM;gi)oo(yMt;g2)^ 

This matrix integral have the perturbation series, described by the following discrete model: 

^ (■?! ;9l)h-_„ + l{92 ;92)h/_„ + i 

Z„(l^;t,t) = (xy) 2 ^0,0 2^ A(gi)A(g2 )lle 

/i,/i'>0 j=l 

For axial-symmetric case (jl.4j) . and for the case q\ = q2 = g, we have the model 



^n(Kxi;t,t') = (x?/)^^?o,o$^(A(g'^))2n^ 

/i>0 i=l 

(4.16) 

where the right hand side in the case |g| = 1 can be interpreted as a discrete one- matrix model 
Z5, and for gi^2 = pe^^^^ it is a version of discrete MNM model Z\. 



(E) Kontsevich-like discrete model Z^. To get model (|2.6|) we consider the axial sym- 
metric case. 

First, we choose the variables t' = t'(?/*^")). Then, we specialize t according to any of ()4.H) . 
(jOI), or jOI). 

Let us consider only the case when t is specialized by ()4.4j) (which is the most general case 
among specializations (jHH), (jOI), (03) )• Then, using (jSUDI), using definition (fTT!^ of 

the Schur function, using ()4.8|) and using antisymmetry of the Vandermonde determinant, we 
can rewrite ()3.3p in the form 



Z„(Kxi,t(a,g),t(yW))A(y) = ^0,0^! ■ J] A(g'^)nk(^Z!!i^l^^ 

(4.17) 

By changing logyi yi we obtain Zq. 

We mark that the relations between continues and discrete models of type, different of 
presented here, were considered in the papers 



5 Dual Schur functions. Associated discrete models and 
discrete matrix models 

In this section, we shall specify the sets of coupling constants t and t' in ()3.3|) in a different 
way, via partitions, say. A* and A'*. To do it, we consider the Schur function S)^{x) in p.l3|) 
as the function of a partition A. We show that if Xi = q^*, where h* are integers, one can 
interchange the roles of x and A. Using this fact we differently relate series (j3.3|l to discrete 
analogs of partition functions of matrix models. 

Given partitions A and A*, whose lengths do not exceed a given number n, we let 

hi = \i — i + n , h* = X* — i + n , (5.1) 

so that we have hi > h2 > hs > ■ ■ ■ > hn > 0, and hi > h2 > h^ > ■ ■ ■ > h*^ > 0. Let 

x*=q''% X, = g'**, i = l,...,n, (5.2) 

where g is a given complex number. 

By definition ()1.13p of the Schur function, we obtain the following Lemma 

Lemma 2 For any q, let x = {xi, . . . , x„) and x* = (x^, . . . , x* ) be as in l[5.^} . Then we have 

A{x)s^{x) = A{x*)ss^,{x*) (5.3) 

Thus, given s^^{x), with x as described above, we obtain the dual Schur function sx*{x*), 
where the roles of h and h* are interchanged. Letting g — 1 in ()5.3|) . we get A{h*)s\(t{n, 1)) = 
A{h)sx*{'t{n, 1)) in accordance with ()4.8p . 

We remark that the case (|4.3|) results from (|5.2p for h* = n — i and n — > 00. 

Given complex qi and q2 and given partitions A* and A'*, we now consider 

Xi = (t , Vi = (I2 , (5-4) 



where 

h* = X*-t + n>0 , h!* = X'*-t + n>0 , (5.5) 

putting 

X'i = Qi , Vi = <l2 > 

= — z + n > , /I'j = A'j — z + n > 

We choose t(x) = (ti, ^2, • • "^'(1/) = (i^'ij '^2; • • •) 

^m = — > gi , ^m = — V?2 , m = l,2,... (5.6) 
m ^ — ^ m ^-^ 

1 = 1 2 = 1 

For the case when Xi = i = 1, . . . ,n, that is A* = 0, we shall write 

1 " 

t(n,g) = (ti(n,g),t2(n,g),...) , t,n(n, g) = - V , m = l,2,... (5.7) 

m ^-^ 

i=l 

Then, using Lemma 1, we obtain for ()3.3|) 

Zr.iV;t{x),t'{y))A{x)A{y) = n\J2det (te) sx4qi)sx'4Q2WQi)^i<l2) , 

\ •' / i,j=l,...,n 

h,h' 

where the sums range over the cones 

oo > hi > ■ ■ ■ > hn > , oo > /I'l >■■■>/;.'„> , 

For the choice of t(x) and t'{y), given by partitions A* and A'* via ()5.6p and ()5.4p . we shall 
use the notation 

Z„(V;A*,A'*;gi,g2) := Z„(l^; t(x), t'(|/))A(x)A(y) , (5.8) 

where the right hand side is defined by ()3.3|) . The case A* = A'* = is related to ()5.7|) . 

Let us note that, in the case n ^ oo, the choice of the coupling constants in form (15.61) is 
rather general. 

Using skew-symmetry properties of the factors, and changing the sum over partitions to the 
sum over all non-negative integers, we get 

1 " 
Z„(y;A*,A'*;gi,g2)=n!<7o,o Yl 7;^A(g^)A(g2'').,.(gf)s,4gr) H (^'^^'^O ' 

h,h'>0 ^ '' i=l 

where we use notation ()5.8|) . and where the sums range over all the n-tuples of non-negative 
integers hi,...,hn and h[,...,h'^. (We kept l/(n!)^ only to explain that the (ra!)^ in the 
denominator is obtained when we change the sum over cones to the sum over all non-negative 
integers) . 

Associated discrete model Let us introduce the following partition function 

Zf-(F;t,t';gi,g2) = ^ Yl A(gf )A(g2'') J] e^^*''''')^^^*'' (5.9) 

■ h,h'>0 i=l 

This sum can be viewed as a partition function of different models, which are discrete analogue 
of the models of random matrices, the choice of the model of random matrices depending 
on the specialization of the complex numbers qi and q2- For instance, for qi = pe^^^ and 
q2 = pe"^^^, where p 7^ 1, sum ()5.9p can be considered as a discrete version of ()1.3p . 

Considering the perturbation series in coupling constants t, t' of (j5.9p we can repeat the 
same trick we used to get Theorem 1. We obtain the series of type (13. 3p . 



Theorem 2 Let 

n 

~9km = AY, n^i''^^"'''"^^"'"' (^^) ' k,m = 0,1,2,... , (5.10) 

h,h'>0 i=l 

where the sum range over all non-negative numbers hi, . . . ,hn and h[, . . . , h'^, and qi and q2 are 
given complex numbers. The perturbation series in the coupling constants t and t' for partition 
function of the model 45. can be presented as follows: 

Zf"(V^;t,t';gi,g2) = h,y{n)s,(t)sy(t') , 

A,A'eP 
e{\),e{\')<n 

where the sums range over all non-negative numbers hi, . . . ,hn and h'l, . . . , /i^, and where 

hi = Xi — i-\-n>0, h[ = X'i — i -\- n > 
This shows that sum ( 15. !J\} is a TL tau function, where n,t,t' play the role of higher times. 

Remark 6 The sums li4-iUj) , H5.9\) can be studied by the method of ( discrete ) orthogonal poly- 
nomials (about the discrete version see, for instance, which links these sums with tau 
functions of integrable hierarchies, as done in TTd^ for the continuous case (namely, when we 
have integrals, instead of sums). 

Consider the differential operator s\*{d), which is evaluated as the Schur function ()1.12j) . 
where each variable tm is replaced by the differential operator dm = {l/'m){d/dim), and where 
the partition A* is related to the variables x by (jS.Sj) and ()5.4|) . Similarly, consider sx'*{d'), 
which is the Schur function ()1.12|) . where each t^ is replaced by d'^ = {l/m){d/di'^), and where 
the partition A'* is related to the variables y by ()5.5|1 and by ()5.4|1 . 

By Lemma 2, we obtain 

Theorem 3 Let g^m = e^*"". Suppose conditions A3, jj)) and / I5.i^) of Theorems 1 and 2 are 
satisfied. Let an '■= go,o{n)/gofi{n) . In the notation 45. we have 

Z„(V;0,0;gi,g2) =a„Zf^'-(V^;t,t';gi,g2)|t,t'=o , (5.11) 

Zr,{V; A*, A'*; qi, q^) = a„Sv id)sx'^ 0') ■ "(V; t, t'; qi, q,) (5.12) 

//, moreover, the potential V is deformed to Vn, subject to the diagram \3. t^) . then we can 
remove the restriction tot, t' = in \5.11\) and h5.12^) : 

Zn{Viy,0,0-qi,q2) = anZf''\V-l,l'-qi,q2) , (5.13) 

ZniYiy, X*, X'*- qi, q2) = ansx4d)sy4d') ■ Z^''''{V; t, V; qi, q^) (5.14) 

We hope that the problem of averaging of the partition function ()1.3|) over initial data, namely, 
over t and t', can be related to the problem of random partitions j32j. A* and A'*. 



Axial-symmetric case (|1.4|) In this case h = h'. We denote V{hi) := Vh^hi- We have 

_ n 

Zf-(y;t,t';g„g,) = ^ A(g^) Afe^) J] e«^*'-^'')+^(^'-^'%^('^^) , (5.15) 



h>0 1=1 

where the sum ranges over all non- negative integer hi, ... , hn, and 



(5.16) 



where t(x) and t'{y) are described by the partitions A* and A'* as described by ()5.6p and ()5.5|) . 
For A* = A'* = 0, we have sx*{d) = sy*{d') = 1, therefore, the partition functions of the 
continues and discrete models coincide. 

We have different models for different qi and q2'. 

(A) If qi = q2 = e'^^^ and all Vmm = 0, then, ()5.15j) . where p = p' = 0, is a discrete version 
of the partition function of unitary one- matrix model J7], [21], that is 



f ^.UTrU^^A^u-^dJJ= J2 ^A(t).A(t') = r.(n,t,t') 

"^^(") A:£(A)<n 



(5.17) 



where we set t = t, t' = t'. In ()5.17|) . the r is the following step function: r{k) = 1 if A; > 0, 
r{k) = if A; < 0. 

(B) Now, qi = q = pe^^^ is a complex conjugate to q2, and |p|, \h\ <^ 1. Let each t'^ = tm 
is a complex conjugate to tm, that is, by (j5.6|) . we have h* = h'* for each i. We obtain 



Z„(Kxi;0,0;g,g) = ^ | A(g'^) f J] ^""^'^^ 



h>0 



i=l 



which is a discrete analogue of the axial-symmetric model of the normal matrices ()1.3|) and 
(ITD . and 



ZniV^,, A*, A-; gi, g^) = ^.a*(5)s,.(9') ■ U e^^''^''^'-^^'''^''^e^^''\i,=o 



h>0 



i=l 



6 Discrete models as soliton solutions 

Here we consider the axial-symmetric case. In Frobenius notations (see ^) , X = {a\(3) , we have 
and the series ()3.10p may be restated as 



^n(Kxi;t,t') = n!5fo,o 



oo oo 



1+E E 



^S(a|/3)(t)s(„|/3)(t') 



V 



k=l ai>->a^.>0 
n>/3i>->/3fe>0 



(6.1) 



/ 



According to ([IH]), when we specialize t by ()4.1|) - ()4.4|) . we obtain soliton tau functions of 
certain (dual) TL equation. 

For instance, if we take t = t^, we obtain 



(^!5'o,o) 005 ^ ) 



k=l ai>-->ak>0 
n>/3i>->/3fc>0 

Then, if we take the collection ^o, (,2, ■ ■ ■ by C,m '■= ^(t, m + a) —^(t' , {m + a)~"^) = + 
(ij'^ik — + 0')~^i'k)y which, for any choice of t', is a degenerate soliton tau function of a (dual) 
TL hierarchy, where the variables tk,t'k play the role of higher times, see for details. 

Another example: if we take t = t(oo,g) and choose the collection ,^0) ^1, ^2, ■ ■ ■ by = 
^(t, g™") — ^(t', q~"^), then we obtain another soliton tau function of the dual TL equation: 

n>/3i>->/3fc>0 

(6.2) 

The partition functions of our discrete matrix models ()4.12|) . ()4.14|) . ()4.15|) . and ()4.16|) are 
also soliton tau functions. This fact is clear after re- writing sums using the Frobenius notations. 

It will be interesting to compare these results with [2H1, where some models of random 
matrices where interpreted as a limit of infinite-soliton tau function in a different way. 

7 Two-matrix models with generalized interaction term 

(A) Development of determinants in the Schur functions Consider the power series 

00 

/(^) = ^e«"-«"z" (7.1) 

n=0 



Lemma 3 Let x = (xi, . . . , n), y = {yi, . . . , ?/„). Then 

•^-tIw^^ ^ (7.2) 

^ ^ AeP;£{A)<n 

where A(x) = ni<i<j<n(^« ^i)' = 1, and c„ = YY^Zq e^o"^" for n> I. 

Putting a„ := e^"~^" so that e^'^i*-^'"*"^""*-* = CnYYi^iCLhii defining n x 00 matrix M by 
Mij = ttjxj and 00 x n matrix by Njk = yl, and using p.l3|) . we observe that ()7.2|) reduces 
to the following well-known formula for the determinant of product of two matrices: 

det(MAr) = ^ det(Mt;;-'^) det(Ar,i-;;:;,J . 

hi>h2> — >h„>l 

While this may be the simplest proof of the lemma, the following proof based on the known 
facts about tau functions is equally simple: First, we identify f{xiyj) with the tau function 

Tr{l, Xi, y,) = 1 + r{l)x,yj + r{l)r{2)x^y] + ■■■ (7.3) 

by r{k) = e^'=~'^'=~i. It may be shown by Wick's theorem (or the determinantal Fay identity |^ 
Theorem 4.2]) that 

det(rr(l, Xi, 'u,))fi._i 
r.(n,t(x),t*(y)) = ^ ^(;);|;) -'-^ , (7.4) 



where 

n-l 



ci=i, c„=n^(^)'""' c^-s) 



fc=i 



It means that for any power series / as in ()7.ip we have 

(Here c~^A{x)A{y) is assumed to be non- vanishing; otherwise formula ()7.4p should be modi- 
fied.) 

Remark 7 Choosing r{k) = {k + n — M)'^ Y[^=ii^ + ^i) Y[i=ii^ + h)^^ for the tau function 
of hypergeometric type, we obtain the hypergeometric functions of two matrix arguments X, Y 



s 



X,F) =r,(M,t(x(")),t*(y("))) 



ai + M, . . . , ap + M 
61 + M, ■ ■ ■ , 6, + M 
^ nLi^A(t(afc + M, 1)) ^ £A(X)£A(r) 

£(A)<n 

i/ere x'-"-' = (xi, . . . , a;„) and y^"-* = (yi, . . . , ?/„) are the eigenvalues of the matrices X and Y , 
respectively. 



Examples of ■ o^o(^, 4) =e^'^ , i^o(a|^, Q = det(J„ - X)~" . Example of (fT^ : 



ai + n, . . . , ttp + 
61 + n, ■ ■ ■ , 65 + n 



c« , / p / ai + 1, . . . , ttp + 1 
A(x)A(y) ^^"^^ + + l 



i,j=l 

(7.7) 

which relates hypergeometric function of matrix arguments to hypergeometric functions of one 
variable. 

(B) Angle integration of determinants We shall exploit the following formulae for inte- 
grations of Schur functions over the unitary group [H]. Let d^U be the normalized Haar measure 
on U{n), and let 6^^x be the Kronecker symbol. By we shall denote by n unit matrix. Let 
A and be n by matrices with the eigenvalues ai, . . . , a„ and 61, . . . , 6„ respectively. (We 
are interested mainly in hermitian and unitary matrices, where eigenvalues are, respectively, 
real numbers or numbers on the unit circle). Then 

s,{AUBU^)dM = , (7.8) 

U{n) Sx{^n) 

and then consider the following model of two hermitian (or two unitary) random matrices 

Jn = j e^'^'^^''^+^''''^^''^g{M,M2)dM,dM2 , (7.9) 

where the measures dMi and dM2 are defined in a standard way (see jSSl); and where the 
interaction term g is a hypergeometric tau function (see p.l7|) ) of the form 

g{X) = r,(n,t(X),t(/„)) := ^ e«'^«-«— Sa(X)sa(/„) , (7.10) 

ASP 

^(A)<n 



where r{k) = e^'= is not fixed, and where 

t(X) = (i Tr X, i Tr X', ^ Tr X^ . . .) (7. 11) 

Then due to ()7.8I1 it is possible to perform the angle integration over f/(n), where Mi = 
UiAUl M2 = U2BUI, A = diag(ai, . . . , a„), E = diag(6i, . . . , 6„). 

/ g{M,M2)dM = Trin,t{Mi),t{M2)) (7.12) 

Ju{n) 

Finely we obtain the following integral over eigenvalues a, and bf. 

Jn = c„/A(a)^A(6)^^^4x#n^^^^"^"'''^^'^^^«^^^^ 
J A (a) A (6) AJ. 

/n 
A(a)A(6) H fiaA)e^'^'''^+'''^''^da,dk (7.13) 
i=l 

8 Conclusion 

Starting with the partition function of the model of normal matrices and using the duality 
relations for Schur functions, we get a Zoo of discrete models. In axial symmetric case these 
discrete models may be identified with soliton tau functions. The next step is to study the 
n — i> 00 limit of these models. Our results means that we should obtain a behavior common 
for groups of different models. Let us notice that discrete models Zi — Z^, we have obtained, 
can be studied by the method of orthogonal polynomials. It will be interesting to consider 
continuous vs. discrete models in the framework of different duality problems, see [31]. These 
topis together with a generalization of all results to multi-matrix models we shall study in a 
more detailed paper. 
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